Chern-number spin Hamiltonians for magnetic nano-clusters by DFT methods 
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Combining field-theoretical methods and ab-initio calculations, we construct an effective Hamil- 
tonian with a single giant-spin degree ol freedom, capable ol the describing the low-energy spin 
dynamics ol lerromagnetic metal nanoclusters consisting ol up to a lew tens ol atoms. In our pro- 
cedure, the magnetic moment direction ol the Kohn-Sham SDFT wave-lunction is constrained by 
means ol a penalty functional, allowing us to explore the entire parameter space of directions, and to 
extract the magnetic anisotropy energy and Berry curvature functionals. The average of the Berry 
curvature over all magnetization directions is a Chern number - a topological invariant that can 
only take on values equal to multiples of one half, representing the dimension of the Hilbert space of 
the effective spin system. The spin Hamiltonian is obtained by quantizing the classical anisotropy 
energy functional, after performing a change of variables to a constant Berry curvature space. The 
purpose of this article is to examine the impact of the topological effect from the Berry curvature on 
the low-energy total-spin-system dynamics. To this end, we study small transition metal clusters: 
Co n (n = 2, 5), Rh 2 , Ni 2 , Pd 2 , Mn x N y , Co 3 Fe 2 . 



I. INTRODUCTION 

The present interest in magnetic nanoparticles, nano- 
clusters and molecular magnets^ 2 - is partly motivated by 
the quest for nano-scale information storage devices^, 
and by possible applications in nano-spintronics and 
quantum computation devices. Most of the advances 
achieved in the last few years in this field have been fu- 
eled by the remarkable experimental ability to synthesize, 
manipulate and characterize these systems at the atomic 
scaled. Theoretical work, based on refined models^^ 
and advanced computational tools^, is also provid- 
ing crucial understanding and new predictions. Mag- 
netic metal nano-clusters, like other magnetic nanostruc- 
tures, typically display enhanced magnetic properties^ 
and many novel classical and quantum phenomena, which 
challenge our understanding of magnetis m 10 ! 11 . 

A quantity that is the topic of intense interest in nano- 
magnetism, is the magnetic anisotropy. This is the de- 
pendence of the total energy of the system on the orienta- 
tion direction of the magnetic moment. The presence of 
two or more stable orientations in the anisotropy energy 
functional, separated by an energy barrier, is the crucial 
property exploited in information storage devices. The 
origin of Magnetocrystalline Anisotropy Energy (MAE) 
in magnetic materials is subtle, being a relatively small 
effect arising primarily from the spin-orbit interaction. It 
is now possible to investigate experimentally the MAE in 
individual magnetic nanoparticles^ and nano-clusters^ 
with unprecedented accuracy. Typically the MAE in 
these magnetic nanostructures is found to be larger than 
in bulk. Gambardella et al. 1A found that monoatomic 
Co chains on a stepped Pt(lll) surface are ferromag- 
netic and have a MAE that is one-order of magnitude 
larger than bulk Co. Understanding the reasons of this 
enhancement and how to increase it further is one of the 



goals of ongoing research. Recent theoretical work pre- 
dicts that the MAE in unsupported 4d transition metal 
(TM) chains^ and dimersi& have giant MAE up to 50 
meV. 

Another topic of great importance, connected to the 
MAE and equally intensively investigated, is the spin 
dynamics of the collective magnetization orientation de- 
gree of freedom in single-domain nanoparticles. For 
a nanoparticle containing several thousand atoms, the 
magnetic orientation degree of freedom is a classical vari- 
able that oscillates around one energy minimum of the 
MAE or switches between two of them by thermal ac- 
tivation. Thermal switching in individual superparam- 
agnetic nanoparticle has been recently studied by spin- 
polarized scanning tunneling microscope^. When the 
nanocluster is sufficiently small, the quantum behavior 
of the collective spin becomes important. The crossover 
from classical to quantum behavior for the magnetiza- 
tion collective variable is a complex problem of great sci- 
entific interest. The physical systems on which research 
has focused so far in studying this question are molecular 
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These systems contain about 



100 atoms per particle, are insulators and behave like a 
single microscopic quantum spin, e.g. S — 10 for Mni2 
acetate. Quantum effects are clearly discernible in their 
low-temperature relaxation properties. A large amount 
of literature has been devoted to studying these systems 2 . 
Similar studies of the quantum spin dynamics in engi- 
neered magnetic metal nanoclusters have just begun. In 
a recent key experiment, Hirjibehedin et al2^ investi- 
gated collective spin excitations in individual scanning 
tunnel microscope (STM)-engineered linear chains of 1 to 
10 manganese atoms, by means of inelastic electron tun- 
neling spectroscopy. Spin excitations that change both 
the total collective spin and the spin orientation were 
observed. 
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On the theoretical front, the microscopic description 
of quantum spin excitations in magnetic metal clusters 
is more challenging than in insulating molecular mag- 
nets, since in the metal clusters these excitations are in- 
tertwined with electronic particle-hole excitations. The 
need to include the crucial spin-orbit interactions adds to 
the complexity of the task. However, in a paper published 
a few years ago^i, two of us put forward the idea that, 
for ultra-small metal clusters, where the single-particle 
energy level spacing is typically larger than then total 
MAE barrier, the "fast" (high-energy) electronic degrees 
of freedom can be integrated out. It is then still possi- 
ble to describe the "slow" (low-energy) dynamics of the 
collective magnetization orientation in terms of a single 
quantum mechanical giant spin J only, like in molecu- 
lar magnets. When spin-orbit interactions are neglected, 
the ground state of a magnetic nano-cluster will tend 
to have a large total spin quantum number S, and may 
have additional orbital degeneracies g or i, as well. In the 
atomic limit for example, the ground state is character- 
ized by total spin and total orbital angular momentum 
L quantum numbers so that g or b = 2L + 1. In clusters, 
orbital symmetries are reduced and g or b will tend to be 
small - in the simplest cases equal to 1. In this case, we 
identify the giant spin J with S. When spin-orbit inter- 
actions are included, this large ground state degeneracy 
will be lifted. Following Ref. [2l|, the approach we take in 
this paper assumes that the magnetic nano-cluster has a 
well-defined group of low- lying energy levels which can be 
identified with the magnetic moment orientation degree- 
of-freedom. For a nano-cluster with total giant spin J, 
2J+1 many-body states will lie in this group. For atoms, 
for example, Hunds third rule states that J = L ± S. 
The spectrum of these levels is described by a giant spin 
low-energy effective Hamiltonian which is the quantum 
generalization of the magnetic anisotropy energy. 

When integrated out, the fast electronic degrees of free- 
dom produce a Berry phase term in the effective action 
that can profoundly modify the behavior of the magne- 
tization orientation degree of freedom. Indeed in this 
approach, orbital contributions, arising from spin-orbit 
interactions, are automatically included via the Berry 
phase term. The giant-spin of the formalism is a topo- 
logical half-integer invariant known as Berry phase Chern 
number, which characterizes the topologically non-trivial 
dependence of the nano-cluster's many-electron wave 
functions on magnetization orientation. 

The purpose of this paper is to carry out the procedure 
suggested in Ref. by calculating Berry phase Chern 
numbers and extracting effective spin Hamiltonians for a 
variety of realistic magnetic nanoclusters containing dif- 
ferent transition metals. In order to carry out this pro- 
cedure, we calculate the MAE for all these systems. We 
treat the many-body electronic problem using state-of- 
the-art numerical methods based on Spin-Density Func- 
tional Theory (SDFT). For transition metals having the 
largest MAE such as Cobalt, our procedure is in prin- 
ciple valid for clusters containing up to hundred atoms. 



However, with current state-of-the-art plane wave SDFT 
codes and the advantages of parallel processing, the com- 
putational load sets our estimated upper limit at a few 
tens of atoms. Here we investigate unsupported ultra- 
small clusters with no more than 5 TM atoms. We 
neglect non-collinear spin configurations and limit our- 
selves to coherent magnetization solution, although we 
allow the possibility of both ferromagnetic and antiferro- 
magnetic order parameters. We focus, in particular, on 
the challenging situation, quite ubiquitous in metal clus- 
ters, where the HOMO-LUMO gap vanishes. This is the 
case where the Berry phase contribution is crucial. By 
comparing our treatment with the ordinary procedure of 
extracting spin Hamiltonians from the MAE functional, 
where the Berry phase term is neglected, we show that 
the latter can lead to incorrect results for the thickness 
of the magnetic anisotropy energy and for the oscilla- 
tion frequencies of the magnetization orientation degree 
of freedom. 

The paper is organized as follows. In section II we give 
a short description of the theoretical background and the 
technical details needed to understand the output of the 
model. We then present our main results in section III 
and proceed with a more in-depth analysis of the dif- 
ferent systems in section IV. Finally, we summarize our 
conclusions in section V. 



II. THEORY 

A. The Berry curvature and the Chern number 

The fundamental origin of the Berry phased is some 
parametrization of a given term in the Hamiltonian, rep- 
resenting a coupling to the external world 54 It is this 
parametrization that leads to an observable which can- 
not be cast into the form of an Hermitian operator, but 
is instead given in terms of a phase change in the sys- 
tems wavefunction. In the systems we are interested in, 
this parameter is the orientation of the coherent mag- 
netization direction of the magnetic cluster, represented 
by the unit vector n. The dependence of the many-body 
wavefunction on the unit vector h gives rise to the Berry 
phase. In its most fundamental form, the phase change 
in the ground state wavefunction as we vary the magne- 
tization between points 1 and 2 is given by 

Aip 12 = -Imlog(*(ni)|*(n 2 )) = -Imlog(l|2), (1) 

which is an arbitrary phase that can be gauged 
away. By closing the loop in traversing the vertices 
of a triangle we obtain a non-vanishing phase change: 
— Im log(l|2)(2|3)(3|l). In the continuum limit of a dis- 
crete, closed path 7 we obtain the usual definition of the 
Berry phase 

V = i<fdh- <*|V ft *). (2) 
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If we then convert the line integral to a surface integral 
over the enclosed area, we can identify the integrand as a 
very convenient, gauge invariant quantity and define the 
Berry curvatur o 23 ' 24 , 



C[n] = iV fi x (tf |V fl tf) 



(3) 



The dimensionality of the total spin system is encoded 
in the dependence of the many-body wavefunction on the 
magnetization direction as a topological invariant called 
the Chern number—, obtained by taking the average of 
the Berry curvature over the unit sphere^, 



"-7- 
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C[n] ■ h dA. 



(4) 



Below we will refer to C[n] = C [n] ■ n as the Berry curva- 
ture. Without SOI, the C[n] is constant and is trivially 
equal to the total spin quantum number S=a. There is 
no anisotropy as the magnetization direction is varied. 
Anisotropy comes from the coupling of electronic spin 
and orbital degrees of freedom via the inclusion of SOI. 
The dimension of the Hilbert space where the systems 
total effective spin resides is now the good quantum num- 
ber and total angular momentum J, replacing S as the 
Chern number of the system. For example, consider the 
case of a single Co atom: disregarding SOI, the system 
has a total spin of S = § coming from the 3d electrons. 
Including SOI, we instead find a Chern number equal to 
J = |, in accordance with the ground state predicted 
by Hund's rules. Of course there is no anisotropy for a 
single atom as it is rotationally invariant in the absence 
of a molecular or crystal field. 

When the size and geometrical complexity of the clus- 
ter increases, determination of the total J quantum num- 
ber quickly becomes a highly non-trivial task. As is com- 
monly the case, one may assume that SOI does not cause 
a level-crossing at the HOMO level. The dimensionality 
of the system can then be inferred from the total mag- 
netic moment in the absence of SOI, after which the spin- 
orbit shifts can be estimated perturbationally. Adding 
the caveat of a dense level structure at the HOMO level, 
as in small transition metal clusters, there is the possibil- 
ity that a spin-orbit induced level-crossing may alter the 
dimensionality, rendering the aforementioned reasoning 
specious. In the presence of SOI, one may again attempt 
to infer the total spin dimension from the magnetic mo- 
ment, but this is a risky process as the numerical diffi- 
culties associated with its determination for TM clusters 
can easily cause a round-off error. By contrast, the Chern 
number can only take on values of multiples of one-half, 
resolving any ambiguity in obtaining the dimension of 
the total spin system. 

It is useful to consider the perturbational expression 
for the Berry curvature. In the so-called parallel trans- 
port gauged 3 - the phase of the wavefunction is kept con- 
stant as the magnetization direction is infinitesimally var- 
ied. Each infinitesimally translated state then becomes 
orthogonal to the previous state. This means that the 



wavefunction can become multiply valued and obtain a 
sign change as h executes a closed loop on the unit sphere. 
The choice of parallel transport gauge is implicitly made 
in the perturbational expression, 



orii 



E 



<* m |^?l*°> 



dru E (n) - E m (n) ' 



(•5) 



where we denote the ground state by ^° and the sum 
runs over excited states ^> m . (0 can be inserted into the 
expression for the curvature (0 to yield 
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i=l m^O 



I dn k 1 



[E (n) - E m (n)Y 



-Tli £ijk- 



(6) 

Although this equation does no service in calculating 
the curvature (one would have to calculate the excited 
states) , it clearly exhibits the strong dependence of the 
curvature on the HOMO-LUMO gap. 



B. Quantum Hamiltonian Extraction 

We follow the procedure outlined in2i and start from 
an approximate imaginary-time quantum action that is a 
functional of the magnetization direction n of the single 
effective spin degree of freedom, 



S[n] 



On 

(*[n]|V ft *[n])- — +E[n] 

OT 



(7) 



The orientational degree of freedom n is more pre- 
cisely defined as the direction in spin-space of the SDFT 
exchange-correlation effective field^i, or its spatial aver- 
age if low-energy states have non-collinear magnetization. 
In Eq. ([7]) is the Kohn-Sham single-Slater determi- 

nant ground state defined by this orientation, and E[h] is 
the SDFT energy. Here spin-orbit terms are included in 
the Hamiltonian. For simple model Hamiltonians where 
the exchange interaction can be decoupled by means of 
auxiliary field functional integrals, Eq. {7} can be derived 
microscopically, starting with a path integral over elec- 
tronic coherent states. Eq. ([7]) was originally proposed 
by Niu and Kleinma n 28 ' 29 to study the spin-wave dynam- 
ics of bulk TM ferromagnets in the context of the adi- 
abatic approximation, and successfully exploited in the 
actual derivation of magnon dispersion by means of den- 
sity functional theor y 30 ' 31 . 

The Euler-Lagrange equations for extremizing the ac- 
tion Eq. (0 yield the semiclassical dynamics of the ori- 
entation direction h, 

dE 
dn a 



-iC, 



'a/3 [n] hp + 







(8) 



where the Berry curvature matrix elements C a p are de- 
fined as 

C aP [n] = ^<*[n]|^-|*[n]) - ^<*[fi]| J?-|*[ft]) . 
on„ onp orip on a 

(9) 
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These equations reduce to the classical Landau-Lifshitz 
(LL) equations when the curvature is constant and equal 
to the total spin moment S. In the presence of SOI 
however, the curvature is in general non-constant and 
can deviate considerably from S, particularly along di- 
rections in correspondence of which the HOMO-LUMO 
energy gap is small. Thus the LL equations can fail in 
this case. In this paper, we aimed instead at deriving 
a quantum mechanical description of the effective total 
spin, whose magnitude we identify with the Chern num- 
ber S. Equipped with the magnetic anisotropy energy 
and the Berry curvature functionals on the unit sphere, 
we can proceed with the extraction of an effective quan- 
tum Hamiltonian for this spin. Assuming that the cur- 
vature is positive definite, the low energy effective action 
(|7|) can be mapped onto that of an effective quantum 
spin Hamiltonian by first performing a variable transfor- 
mation to constant curvature space. This Hamiltonian is 
a convenient tool that can be used for calculations of tun- 
neling amplitudes, non-linear response to external fields 
etc. The variable transformation is given hy^L 



(10) 



u — > u = — I 
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j du"J d<t>"C{u", <£")• (ii) 



with respect to the spin coherent states |s, n' (i')} describ- 
ing a spin s parametrized by the unit vector n' (t'). 

The quantization of the transformed magnetic 
anisotropy energy is then performed by expanding in 
spherical harmonics of highest possible order 2s, where 
time-reversal symmetry precludes odd-powered harmon- 
ics from contributing. 

E[h(h'(t'))} = f^ E TYl n {n') (15) 

1=0 m=-l 

The spin coherent state \s,n' (t 1 )) can be constructed by 
applying rotation operators in the single spin-space act- 
ing on the coherent spin state with maximum polariza- 
tion m — s. This allows us to express the RH side of (TT4")) 
in terms of spherical harmonics. 

E {«, ^ (R) k m ) ( s , m\H\s, m!) (s, m'\V (R) \s, s) 

mm' 

= E ( R lm) Kmm'Rm's ( i6 ) 
mm' 

= E (-i) m ~VM2Z + i) 

(17) 



mm'lrh 

s s I \ / s s I 



\ I n I I / - I Yi % (n) H s mm 

-s s v I \ —m mm 1 1 v ' mm 



where C is the curvature, iS the Chern number and u = 
cos 9. This change of variables rescales the local curvature 
metric such that 



C (u,(j>) dud<j) = Sdu't 



(12) 



From this expression we see that the transformed in- 
finitesimal area element will be rescaled in order to com- 
pensate for variations in the curvature. Consequently, if 
C > S the variable transform will cause an expansion of 
the grid points and if C < S a contraction. At each grid 
point there is an energy value associated with the mag- 
netization direction that will be relocated by the trans- 
form. In particular, we see from expression ^ that if the 
HOMO and the LUMO level become quasi-degenerate for 
one or more magnetization directions, the curvature will 
become quasi-singular, causing a substantial deformation 
of the grid. 

The variable transformation allows us to rewrite the 
real-time action as 



^spin \J 



dt' 



iA.^-E[h(h'V))] 



(13) 



where A — scf>' (I — cos (9') /sin#'. This is the quantum 
action for a total spin of quantum number s, with a clas- 
sical Hamiltonian H [n 1 (t 1 )} = E [n (h' (t'))^, that is 
equal to the expectation value of a quantum spin Hamil- 
tonian H, 



H[n' (t')} = ( S ,n' (t) \H\s,n> (t')>, 



(14) 



where we have made use of the Wigner-Eckhart theorem 
and expressed the Clebsch-Gordan coefficients in terms 
of Wigner-3J symbols. Equating (TIT)) to the RH side of 
(|15p . convoluting with and inverting, we obtain 



H s — 



(-ir" s E^ 

A fx 



2A + I 



s 
rn 



■in 

s 
s 



s X 

-s 



(18) 



Equation (jTSJ) is of great practical use as gives an an- 
alytical relation between the spherical harmonic expan- 
sion coefficients and the matrix elements of the quantum 
Hamiltonian. Once the Hamiltonian matrix has been ob- 
tained, it can be decomposed in terms of the spin opera- 
tors S x , S y , S z and S 2 . 



C. Computational Details 

Our numerical tool is SDFT using the Projector Aug- 
mented Waves (PAW) formalism 32 as implemented in the 
VASP— code. As such, the method depends on repre- 
senting the atoms with PAW pseudopotentials, in which 
the core electrons are described on an auxiliary radial 
grid. This method can be said to be of all electron type, 
although only the valence electrons participate in the full 
variational procedure, whereas the core region is adjusted 
by applying boundary and conservation conditions. The 
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pseudopotentials are generated in the generalized gradi- 
ent approximation which gives a better description of the 
exchange-correlation effects. 

We use accurate parameter settings for our calcu- 
lations. By this we mean a valence (auxiliary) en- 
ergy cutoff of 348(621), 298(488), 351(709), 326(541), 
381(751), 351(740), 520(815) eV in the case of Co, Rh, 
Ni, Pd, Fe, Mn and N, respectively. In that order, the 
PAW pseudopotentials 3 ^ for these elements are gener- 
ated in the following electronic valence configurations: 
3d 8 4s 1 , 4d 8 5s\ 3d 9 4s\ 4d 9 5s\ 3p 6 3d 7 4s 1 , 3p 6 3d 6 4s 1 
and 2s 2 2p 3 . To avoid wrap-around errors we chose a 
dense FFT-grid taking into account all G- vectors that are 
twice as large as the vectors in the basis set for the valence 
region and the 16/3 times as large for the auxiliary re- 
gion. Interaction between the periodically repeated cells 
is suppressed by selecting a relatively large supercell-size 
of 11 x 11 x 11 A (in the case of Mn x N y we use elongated 
cells suitable for these geometries). Only the T-point is 
used in the calculations and we have monitored disper- 
sion by comparing the eigenenergy spectrum at the T- 
point to that of the X-point, at the edge of the Brillouin 
zone - in general there is no or very little dispersion. By 
making use of the Vosko-Wilk-Nusair— interpolation for- 
mula for the correlation part of the exchange-correlation 
functional, magnetic moments and energies are enhanced. 
All calculations are performed in non-collinear mode (as 
described in ref. [H) with and without SOI. 

The calculation of electronic properties for small TM 
clusters is a daunting task. Well converged solutions 
are often obscured by a dense level structure at the 
HOMO level, causing sporadic jumps between different 
spin-configurations and in a faulty ordering of the eigen- 
levels. As noted by Pederson et. Al21 (as well a o 38 ' 39 ), 
the use of a finite (in our case Gaussian) smearing makes 
it possible to achieve convergence in the self-consistent 
state. This fake temperature causes a decreasing frac- 
tional occupation of the levels around the HOMO level, 
enabling the correct ordering of eigenstates. At the end 
of the calculation, the limit of zero smearing is taken, 
eliminating the artificial entropy contribution to the to- 
tal energy. 

SOI is included and the pseudopotentials are generated 
in the scalar relativistic approximatio n 4 ^ 41 . In this ap- 
proach the relativistic Dirac equation is separated into 
two parts. The first part has a j-weighted/averaged 
scalar-rclativistic Hamiltonian where the dependence on 
k = ±(j + 1/2) has been removed. The second part has 
a spin-orbit Hamiltonian, 

Hbo = W^^ (19) 

where s and £ are the spin and angular momentum op- 
erators and the relativistically enhanced electron mass is 
given by M — m + (e — V) /2c 2 . Here V is the effective 
potential, c the speed of light and e the eigenenergy. The 




FIG. 1: The vertex skeleton employed. At each vertex the 
ground state energy and the Kohn-Sham wavefunction are 
obtained and for every triangle we may form the three over- 
laps needed to calculate the curvature. 

scalar-relativistic Hamiltonian contains the mass veloc- 
ity and Darwin corrections. It is solved by diagonaliza- 
tion, whereafter the spin-orbit Hamiltonian is included, 
such that the full Hamiltonian is solved using the scalar 
relativistic wave functions as basis set. This approach 
gives spin-orbit shifts that are in excellent agreement 
with experimental values, although for the heavier el- 
ements the p-clcctron j-wcighting/averaging procedure 
introduces an approximation error—. 

Our solution for a given magnetization direction is sta- 
bilized via the introduction of a penalty functional in the 
SDFT Hamiltonian, the additional total energy contribu- 
tion of which can be made vanishingly small. This proce- 
dure allows us to calculate the total ground state energy 
for any given direction on the unit sphere. In addition, 
we obtain the Kohn-Sham wavefunction corresponding to 
the particular magnetization direction, which can be used 
to calculate overlaps between wavefunctions for different 
magnetization directions, and subsequently extract the 
Berry curvature. The phase obtained in moving around 
a closed loop 7 may be written in terms of the occupied 
Kohn-Sham single particle orbitals ipi as^ 3 - 

N 

V = iJ2f( ( Pi(n)\Vh^i(n))dl. (20) 
1=1 7 

In practice, we cover the unit sphere with small trian- 
gles and discretize (|20|) to a closed path with three ver- 
tices, expressing it terms of the determinant of the over- 
lap matrices Oij between the Kohn-Sham wavefunctions 
obtained at each calculation vertex, 

Vi = -Im log (det 0\ 2 det 0\ z det O^) . (21) 

where the index i labels the triangles. Dividing (|2"Tj) by 
the subtended area, we obtain the local curvature values 
Ci. Averaging these over the unit sphere yields the Chern 
number, which can only take on multiples of one-half. 

The computational load of this scheme is considerable 
and it is therefore necessary to choose an optimal cov- 
ering of calculation vertices on the unit sphere. A con- 
venient geometry is the 4-frequency icosahedron. An n- 
frequency icosahedron is an icosahedron in which each 
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FIG. 2: The considered clusters with distances in Angstrom. 
The asterisk indicates non-ground state configurations. In the 
Mn 3 NJ cluster the atoms are equidistantly spaced at 1.8 A, 
like in Mn 3 N 2 . 



of the 20 constituent triangular surfaces has been sub- 
divided into n 2 new triangles (see for example ref. l43l ). 
Fig-fflshows the 4-frequency icosahedron with its 162 ver- 
tices, at which the ground state energies the associated 
Kohn-Sham wavefunctions are calculated, and 320 trian- 
gles, for which the Berry curvature values are extracted. 



III. RESULTS 
A. Structural Optimization 

Structural optimization of the Co n (n = 2, ...,5), N12 
clusters is carried out starting from initial atomic posi- 
tions determined by Castro et. Al4i, whereby the in- 
teratomic forces are minimized and the structures re- 
fined using an accurate quasi-Newton method. Initial 
distances for Pd2, Rti2 were taken from ref. l45l In ac- 
cordance with the substrate-enforced geometry (as in 
ref. |20[ ). the Mn atoms are frozen equidistantly at 3.6 
A, with N in between. The Mn^N* clusters are marked 
with an asterisk to indicate that the positions have not 
been relaxed. We will also examine the effect of allow- 
ing the atoms to relax their position on the axis. For 
Mn2N3 we start relaxation from an initial state where 
the N atoms are in plane but off the Mn axis. 

In addition we consider the equilateral Co trimer 
(marked C03) and the effect of reducing the dimer dis- 
tance in Pd2 (marked Pdj). The considered structures 
are shown in fig. [2] where the distances are given in units 
of Angstrom. Table[T]shows the resulting symmetries and 
binding energies. 



System 


Symmetry 


E B [eV] 


S 


J 


Rh 2 


C 2 h 


3.68 


2 


4 


Ni 2 


C 2 h 


2.87 


1 


1 


Pd 2 


C 2 h 


1.36 


1 


1 


Pd* 2 


C 2 j, 


0.69 


1 


2 


Co 2 


C 2 j, 


3.57 


2 


4 


Co 3 


C 2 h 


6.46 


7/2 


7/2 


Co* 3 


Czh 


6.31 


5/2 


3/2 


C04 


D 2d 


10.01 


4 


4 


Co 5 


D Ah 


14.55 


13/2 


13/2 


Co 3 Fe 2 


D 3h 


14.22 


13/2 


11/2 


Mn 2 N 3 


Clh 


14.66 (14.69) 


5/2 (1/2) 


5/2 (1/2) 


Mn 2 N 3 


C 2 h 


12.39 (12.40) 


5/2 (1/2) 


5/2 (1/2) 


Mn 3 N 2 


C 2 h 


12.83 (13.36) 


13/2 (5/2) 


13/2 (5/2) 


Mn 3 N 2 


C 2 h 


12.67 (13.02) 


15/2 (5/2) 


15/2 (5/2) 


Mn 3 N 4 


C 2 h 


20.23 (20.06) 


9/2 (3/2) 


9/2 (3/2) 


Mn 3 N4 


C 2 h 


19.57 (19.36) 


9/2 (3/2) 


9/2 (3/2) 



TABLE I: The symmetries, binding energies and Chern num- 
bers of the considered clusters. The Chern numbers without 
(S) and with SOI (J) may differ in the presence of an SO- 
induced level-crossing. Values in parenthesis refer to the AFM 
configuration. 



B. Chern Numbers 

The Chern numbers for the considered clusters, as ob- 
tained by taking the average Berry curvature over the 
unit sphere, are displayed in table [U 

Fig. [3] shows the 'Chern spectrum' for a few chosen 
clusters. These diagrams indicate the accumulated Chern 
number (the solid line) , obtained by including the single- 
particle Kohn-Sham orbitals up to the given level in the 
wavefunction overlaps, and the level Chern number (the 
points), obtained by taking the difference between accu- 
mulated Chern number of the given level and that of its 
antecedent. We shall refer to the associated quantities as 
the level and accumulated Chern number and curvature 
respectively. 

In Fig. [3l level Chern numbers are formed by an or- 
bital part and a spin part obtained by taking the aver- 
age curvature over the unit sphere. In the absence of 
SOI, the level Chern numbers can only take on values 
of ±i. As the spin-orbit coupling strength is increased 
from zero, only levels that cross can deviate from this 
value by acquiring an orbital part. In accordance with 
the Chern number sum rule for level-crossings found in 
ref. [2l] - the sum of the level Chern numbers before and 
after the crossing is preserved. Notice that the level 
Chern numbers of the dimers are particularly affected by 
level-crossings, attributable to their inherent high degree 
of symmetry. Comparing the Jahn- Teller distorted C03 
with the equilateral trimer C03, we see that the sym- 
metrization of the cluster leads to multiple deviations 
from ±i. A similar situation occurs in C05 when the two 
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FIG. 3: The Chern spectrum for a few selected clusters. The 
points show the level Chern numbers and the solid line the 
accumulated Chern numbers. The HOMO level and the total 
system Chern number have been marked with a vertical line 
and a circle. Deviations from ±| indicate the presence of a 
level-crossing as a function of increasing spin-orbit coupling 
strength. Multiples of half-integers have been marked with 
horizontal lines in order to lend support to the eye. 



axial Co atoms are substituted for Fe. The level-crossing 
effect may in isolated cases even cause neighboring level 
Chern numbers to take on integer values. 

In fig. 21 we plot individual level Chern numbers to- 
gether with their corresponding maximum and minimum 
curvature values on the unit-sphere as error bars, for 
a few representative clusters. Large fluctuations in the 
level curvature signal the presence of a near degeneracy 
between two adjacent levels. In general we find that a 
level with large fluctuations in the Berry curvature is fol- 
lowed by another level with curvature fluctuations of the 
opposite sign that nearly compensate the previous ones, 
(see for example the case of C02 where the curvature 
landscape of two consecutive levels is shown.) This pat- 
tern implies that although the Berry curvature of individ- 
ual levels might be singular in some directions, their sum 
is in general a smooth function of n due to the cancella- 
tion of the opposite contributions of adjacent levels. An 
exception to this rule occurs when the quasi-degeneracy 
involves the HOMO and LUMO levels. In this case, 
the fluctuations of the HOMO are not compensated by 
the unoccupied LUMO; the total Berry curvature of the 
system is completely dominated by the HOMO curva- 




T 3 5 ^1 5 

<Ei>unit Bphera [eV] 



4 -3 -2" H 5" 
^Ep-unilsphera [SV] 



FIG. 4: (Color online) The Chern number plus the maxi- 
mum and minimum level curvature as a function of the av- 
erage eigenenergy taken over the unit sphere. Note the large 
positive maximum at the HOMO-level (set to zero) of C02 
and the induction of high-curvature variations obtained by 
substituting the axial Co atoms in C05 with Fe. 



ture and can therefore deviate strongly from the total 
Chern number and be singular for particular magnetiza- 
tion directions. If the total curvature becomes infinite, it 
may be difficult to implement the procedure outlined in 
Sec. Ill Bl . In particular if the HOMO level has a predom- 
inantly minority-spin character, the total curvature can 
even become negative along these degeneracy directions, 
and the transformations given in Eqs. [TUl [TT]is invalid. 
Physically this means that when HOMO-LUMO gap goes 
to zero the low-energy quantum dynamics of the system 
cannot be described by an effective giant-spin Hamilto- 
nian only. In this case it necessary to include explicitly 
all the electronic degrees of freedom involved in the de- 
generacy at the Fermi level. We will address this issue 
in more detail in Sec. IIVI Depending on the given situa- 
tion, one might still to a good approximation evaluate the 
topological effect when there are only isolated negative 
dips in the curvature. 

In table HI we see that there can be a change in the 
the Chern number, namely the dimensionality of the to- 
tal spin system when SOI is turned on. In general this 
happens when the HOMO-LUMO gap is zero in the ab- 
sence of the SOI. When SOI is turned on, the degeneracy 
is lifted but in a few nearly avoided level-crossings for 
particular directions of the magnetization. Berry phase 
contributions coming from paths encircling these quasi- 
degeneracies are responsible for the change in the Chern 
number. Thus the assumption that the dimension of ef- 
fective spin the does not change when SOI is turned on 
is potentially a hazardous one in the case of small TM 
clusters. 



C. Magnetic anisotropy energies 

Let us first take a look at the magnetic anisotropy 
energies, before performing the variable transformation 
that takes us to constant curvature space. A qualitative 
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System 


Em a/ atom [meV] 


Type 


Rh 2 


26 


EA 


Ni 2 


3.8 


EP 


Pd 2 


1.2 


EP 


Pd 2 


17 


EA 


Co 2 


14 


EA 


Co 3 


1.3/0.7 


EA 


Co^ 


0.8 


EA 


Co 4 


0.9 


EA 


Cos 


0.1 


EA 


Co3Fe 2 


0.5 


EA 


Mn 2 N 3 


0.9/0.4 (0.8/0.4) 


EA(EA) 


Mn 2 N 3 


0.8/0.2 (0.7/0.1) 


EA(EA) 


Mn 3 N 2 


1.4 (0.8) 


EP(EP) 




1.2 (1.0) 


EP(EA) 


Mn 3 N 4 


0.8 (1.2) 


EP(EP) 


Mn 3 NJ 


0.8 (0.3) 


EP(EA) 



TABLE II: The magnetic anisotropy energies per atom. EA 
signifies an easy axis in parallel to the symmetry axis of the 
cluster, meaning the axis perpendicular to the atomic plane 
for Co 3 , the pyramid axis for C05 and Co 3 Fe 2 , the axis run- 
ning through the mid-point of the two opposite solid lines in 
fig- ED for C04. EP signifies a quasi-easy plane perpendicular 
to the symmetry axis running through the dimers (or the line 
connecting the Mn atoms). Co 3 has a high and a low barrier 
separating the global minima that occur when the magnetiza- 
tion is parallel to the elongated side. The low barrier runs in 
a cross-section perpendicular to the trimer plane, whereas the 
high barrier is in a cross-section in the trimer plane. In Mn 2 N 3 
the global minima are perpendicular to the cluster plane. A 
high barrier is separating the oppositely oriented states in a 
plane parallel to the Mn-Mn axis, and a low barrier separates 
the minima in a plane perpendicular to the Mn-Mn axis. For 
the Mn^Ny clusters, the anisotropy is per Mn atom. Values 
in parenthesis refer to the AFM configuration. 



description of the magnetic anisotropy landscapes and 
the associated anisotropy energies per atom, \E max {h) — 
E m in(n)\/N at oms, are presented in table HH 

To evaluate the topological effect, we first emulate a 
simplified approach to construct a quantized Hamiltonian 
by omitting the variable transformation that takes us to 
constant curvature space. We then take the dimension 
of the effective spin-system to be that inferred by the 
magnetic moment in the absence of SOI (i.e. the Chern 
number without SOI in table |T|), with the assumption 
that turning it on causes no level-crossing. The magnetic 
anisotropy energy is taken to be that which we obtain by 
actually switching on the SOI, but should in this context 
be regarded as originating in a perturbational approach. 
Fitting the anisotropy energy with spherical harmonics 
limited by the system dimensionality without SOI, results 
in the effective quantum Hamiltonians displayed in table 



System 


Hamiltonian 


Coefficients [meV] 


Rh 2 


o[l + PSl + 7 S Z 4 ] S=2 


a = 85.1,/?= -0.630 
7 = 0.095 


Ni 2 


q[1 + /3S 2 ] s =i 


a = -7.68,/?= -2.0 


Pd 2 


a[l + pS 2 z ]s=i 


a = -2.44,/? = -2.0 


Pd* 2 


a[l + /3S 2 ] s =i 


a = 66.92, /? = -1.0 


Co 2 


all + BS 2 + 'YS' 4 lq-2 


a = 45.3,/?= -0.555 
7 = 0.076 


Co 3 


a[l + (3S 2 + 

"f{S 2 - Sy)] S =7/2 


a = 0.491,/? = 0.554 
7 = 0.10 


Co 3 


a[l + PS2] s=s /2 


a = 2.88,/? = -0.16 


C04 


a[l + /3S 2 ]s=4 


a = 3.94, /3 = -0.063 


Co 5 


a[l + BS 2 ] s=13/2 


a = 0.589, 13 = -0.024 


Co 3 Fe 2 


a[l + f3S 2 ] s =i 3 /2 


a = 2.71,/? = -0.024 


Mn 2 N 3 


a[l + (3S 2 + 

1\pa — a y)iS=5/2(l/2) 


a = 0.0256(0.73) 
/? = 10.8(0) 
7 = -2.91(0) 


Mn 2 N 3 


a[i + po z + 

l(S 2 — Sy)] S=i /2(l/2) 


a = -0.157(0.50) 
/? = -1.71(0) 
7 = 0.23(0) 


Mn 3 N 2 


a[l + /9S?]5=i3/2(6/2) 


a = -0.358(-0.619) 
(3 = -0.308(-0.800) 


Mn 3 N 2 


a[l + /9Sf]s=i8/2(6/2) 


a = -0.257(3.82) 
/? = -0.267(-0.160) 


Mn 3 N 4 


a[l + /3Sf] S=9 /2(3/2) 


a = -0.303(-1.82) 
P = -0.444(-1.33) 


Mn 3 N 4 


a[l + /3Sf]s=9/2(3/2) 


a = -0.300(1.28) 
f3 = -0.444(-0.444) 



TABLE III: The Hamiltonians obtained by quantizing in the 
spin Hilbert space of the dimension given by the Chern Num- 
ber without SOI (see table U}. These Hamiltonians are calcu- 
lated without performing the variable transformation to con- 
stant curvature space, thereby neglecting the topological ef- 
fect. This approach is designed to emulate the results one 
would obtain by adding the spin-orbit shifts perturbationally. 
The quantization axis of the big spin is taken to be the cluster 
symmetry axis as defined in table [TTT 



der and in isolated cases (C02, RI12) a small fourth or- 
der contribution in the spin operators. The higher order 
contribution is an indication of level-crossings and that 
non-trivial topological effects are present when SOI is 
switched on. If the dimension of the Hilbert space does 
not change when SOI is turned on, there is most likely no 
level-crossing at the HOMO. With a large HOMO-LUMO 
gap, the curvature is essentially constant with very small 
variations around the Chern number and there is no topo- 
logical effect that enters through the variable transform 
(equations [TU] and [TTj) . In this case our effective Hamil- 
tonian approach will yield the same result as in tablc lTTTl 



Generally, we find only coefficients up to second or- 
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System 


Transformed 
Hamiltonian 


Coefficients [meV] 


Rh 2 


a[i + po z + 7^ z 
+SS 6 Z +eSf] J= 4 


a = 298.79, j3 = -2.3602 
7 = 0.94995,5 = -0.11266 
e = 0.0038916 




a[l + 
+7Sj]j=a 


a = 6.19, p = 12.2 
7 = -3.11 


Co 2 


a[l + PSl + ySt 

4-SS 6 + FS S ] j t 


a = 145.30, ,3 = -2.3058 
7 = 0.93035,5 = -0.11038 
e = 0.0038124 


Co* 3 


a[l + /3S 2 ]j= 3 /2 


a = 3.03,/3= -0.445 


Co 3 Fe 2 


a[l + /?S 2 
+7SS/=ii/ 2 


a = 2.61, /3 = -0.0235 
7 = -0.000317 


Mn 2 N 3 


"I* Jj=l/2 


a = 0.137 


Mn 2 N3 


«[5 2 ]jfl72 


a = 0.163 


M113N2 


a[l + /?S 2 

+7J Z J.J=5/2 


a = 2.76, P = 0.365 
7 = -0.0840 


Mn 3 N4 


a[l + (3S z \ J=3/2 


a = 1.71, /3 = -0.444 



TABLE IV: The Hamiltonians obtained by quantizing after 
performing the variable transformation to constant curvature 
space. The cases in table ITTT1 where the transformation exhib- 
ited a very weak and negligible topological effect have been 
omitted above. Only the AFM Mn^Ny have non-trivial cur- 
vatures. 



D. Effective Chern- number spin Hamiltonians 

We now consider the effect on the magnetic anisotropy 
obtained in performing the variable transformation to 
constant curvature space. Fig. [5] shows the total system 
curvature landscapes on the unit sphere for a few clusters, 
as well as the inverse of the HOMO-LUMO gap squared. 
This figure demonstrates that the quasi-singular behavior 
of the Berry curvature landscape is governed by the gap, 
as expected from the perturbational expression Eq. (|6|). 

Table IIVI displays the result of quantizing the classi- 
cal Hamiltonian in a spin-space of a dimension dictated 
by the Chern number with SOI, using the transformed 
magnetic anisotropy energy, thus incorporating the topo- 
logical effect that enters via the Berry curvature. 

Fig.[S]shows the magnetic anisotropy landscapes on the 
unit sphere before and after the transformation to con- 
stant curvature space for Rb.2 and Mn2N3, as well as the 
Berry curvature that enters their transform. Fig.[7]shows 
cross-sections of the unit sphere with transformed mag- 
netic anisotropy energies and comparative sections with 
the untransformed anisotropies, as well as the Berry cur- 
vature for a few clusters. The negative dips at the equa- 
tor for C03 and Co3Fe2 have been cut in the transform 
in order to approximately capture the topological effect 
in these systems. Note that the transform cannot really 
change the magnitude of the energy values, rather, the 
quantization process involves the best possible fit with 
the allowed spherical harmonics which will average the 




FIG. 5: (Color online) The Berry curvature and the inverse 
of the HOMO-LUMO gap squared on the unit sphere for se- 
lected clusters. To the right of each curvature is the inverse 
of the associated HOMO-LUMO gap squared, demonstrating 
the correspondence with the perturbational expression ||B). 
RI12 and C02 display very similar curvatures that are constant 
at 2, with the exception of an equatorial quasi-singular ridge 
that extends up to approximately 16. The AFM Mn2N3 and 
Mn2N 3 have two Gaussian peaks perpendicular to the cluster 
plane, extending from to around 8. The AFM Mn 3 N2 sticks 
out as the correspondence with the HOMO-LUMO gap is not 
entirely obvious. The three lower clusters show the triggering 
of a quasi-degeneracy between the HOMO and LUMO levels, 
by reducing the interatomic distance in Pd 2 , by symmetrizing 
the cluster geometry C03 , and by replacing the axial Co with 
Fe in Co3Fe 2 . 
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FIG. 6: (Color online) The effect of the variable transform 
to constant curvature space for RI12 and Mn2N3. (a) and 
(d) show the calculated magnetic anisotropy landscape on the 
unit sphere and (c) and (f ) show the transformed anisotropies 
where the topological effect of the Berry curvature (panels 
(b) and (e)) has been included. RI12 and C02 are topologi- 
cal twins and have the same curvatures and anisotropy land- 
scapes (anisotropies differ only in magnitude due to the heav- 
ier Rh mass and stronger SOI). The equatorial quasi-singular 
ridge in RI12 and C02 enters the transformation and causes the 
effective barrier to widen significantly. The extreme Gaussian- 
shaped quasi-singular peaks in the curvature of Mn2N2 and 
M112N2 compress the low barriers, so that the transformed 
anisotropy essentially becomes a shallow quasi-plane. 
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FIG. 7: ^-dependence of the curvature and the calculated and 
transformed magnetic anisotropy landscapes for selected clus- 
ters. The thick solid line is the curvature, the thin solid line 
the untransformed anisotropy energy and the dashed line the 
anisotropy energy after the variable transformation to con- 
stant curvature space has been performed. The horizontal 
solid lines indicate the Chern number without SOI, and the 
dashed ones with SOI. Pd and Ni have the same valence con- 
figuration and (although with different HOMOs) are topolog- 
ical twins. 



transformed anisotropy energy into the subspace whose 
dimension is given by the Chern number. As a result of 
this averaging procedure the maximal energy is affected 
in accordance with the effect of the transform entering 
through the curvature. In the case of C03 and Co3Fe2 the 
pull of the sub-Chern curvature on the equator is the pre- 



vailing effect, yielding a lower effective maximum. Not 
very surprisingly, there is no topological effect present 
when the HOMO-LUMO gap grows as in the case of Pd 2 , 
M2, 0034,5 and the relaxed linear Mn x N y (these have 
therefore been omitted from table HV| ). We will now ad- 
dress each group of systems and the mechanisms present 
in detail. 



E. Dimers 

Transition-metal dimers have an inherent symmetry 
that gives them remarkable magnetic properties - they 
are rotationally invariant around the dimer axis. This 
can cause magnetic anisotropy in transition-metal dimers 
to appear already at first order in a perturbational treat- 
ment of the SOt!£ Depending on the electronic config- 
uration and the situation at the HOMO, this symmetry 
property can cause the magnetic anisotropy to be anoma- 
lously strong. 

We will discuss the dimers in terms of their Kohn- 
Sham orbital energies labeled using standard molecular 
notation— Fig. [8] (a) and (d) show these for Rh 2 and 
C02 (around HOMO only) without SOI for majority (|) 
and minority-spin (|), and with SOI and in the hard (x) 
and easy (z) direction. The s and d valence orbitals on 
one atom hybridizes with orbitals on the other atom that 
have the same azimuthal angular momentum me, to form 
bonding and antibonding combinations. In the absence 
of SOI, these are doubly degenerate and are labeled ac- 
cording to their azimuthal angular momentum, where a 
refers to me = 0, 7r to mi = ±1 and 5 to mi — ±2. A * in- 
dicates an antibonding orbital and u/g labels odd/even 
states. The SDFT exchange potential lowers majority- 
spin orbital energies relative to those of minority-spin 
type. 

Both RI12 and C02 have very large magnetic 
anisotropies of 26 and 14 meV per atom respectively. 
They both exhibit an extreme topological effect due to 
the presence of a quasi-degeneracy at the HOMO level. 
The curvature (see fig. [6|) is constant at 2 (note that the 
Chern number is 2 in the absence of SOI) all over the 
unit sphere, except at the equator where a quasi-singular 
ridge extends up to approximately 16. Fig. [8] (b) and 
(f) show the orbital mixing character as given by the 
site projected spherical harmonics inside non-overlapping 
atomic spheres taken over the spinors in the direction 
of the magnetization, for the topmost eigenlevels in the 
hard and easy direction (as the spheres do not overlap, 
some interatomic charge will be missed by the projec- 
tion spheres). Using the calculated Chern spectrum (see 
fig. [3J in combination with the atomic site-projected or- 
bital character for each quasi-particle level, we can ex- 
tract the orbital and spin contribution parts to the level 
Chern number. The crucial level in both Rli2 and C02 
is the HOMO 5*, which is doubly degenerate but singly 
occupied in the absence of SOI. From the Chern spec- 
trum of C02 (fig. 0) and the orbital mixing (fig. [5]) we see 
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FIG. 8: (Color online) Kohn-Sham orbital energies and pro- 
jections of orbital character for RI12 and C02. The left panel 
in (a) and (d) shows the majority and minority- spin levels in 
the absence of SOI and the right panel the levels in the hard 
(x) and easy (z) direction with SOI. The origin of the large 
anisotropy in RI12 and C02 is a first order contribution in the 
spin-orbit perturbation series, resulting from the the splitting 
of the singly occupied but doubly degenerate HOMO-LUMO 
$*-level. The HOMO (number 18) is connected with a dot- 
ted line and its population is marked with a minority spin 
arrow. When the SOI is turned on, the quasi-degeneracies in 
the hard direction split as we change the magnetization direc- 
tion to the easy (symmetry) axis. Panels (b) and (f) show the 
orbital mixing character as given by the atomic site projected 
spherical harmonics taken over the spinors in the indicated di- 
rection of magnetization. Panels (c) and (e) show the orbital 
mixture landscapes of the HOMO on the unit sphere. The or- 
bital mix of the HOMO (18) and LUMO (19) show that along 
the easy axis these states have good mt quantum numbers ±2, 
allowing for a continuous electron path in a plane perpendicu- 
lar to the magnetization direction and a total lowering of the 
energy. 



that the HOMO (level 18) has a level Chern number of 
3/2 with an orbital and spin contribution of 2 and -1/2 
respectively. The orbital angular momentum is opposing 
the magnetization and this level will therefore be shifted 
down, whereas the LUMO-level with orbital part -2 and 
spin part -1/2 will be shifted up. The shifts in the sub- 



levels will to a large extent cancel each other out, but 
the HOMO shift remains uncompensated, causing a very 
large anisotropy energy. In RI12 the spin-orbit shift is 
so strong that the d* descends below the minority-spin 
4s(cr g ), which is much higher in energy for RI12 than for 
C02 - consistent with the larger dimer separation. 

Up to double-counting corrections, the total energy is 
simply the sum of the spin-orientation-dependent shifts 
in occupied orbital energies 5 . In the hard direction spin- 
orbit shifts are small compared with the scale of the 
d-electron spin-orbit interaction (Hso = S,dS ■ L, with 
^° ~ 85 meV and (,f h ~ 140 meV). With the spins po- 
larized in the hard (x) direction the expectation value of 
L x over the unperturbed states is zero, and any shifts 
originate in higher order perturbation terms (at most 
~ 0.1£d, where Wd is the typical d-orbital bond- 
ing energy). However, when the spins are polarized 
along the easy (z) direction, i.e. the dimer axis, the 
first order contribution is non-vanishing and is given by 
±£,d(L z ) = im^/2. For RI12 and C02 the first or- 
der contribution is larger than the obtained anisotropy 
energy and the expected |cos(0)| appearance is not ob- 
served. This is because the cusp at the degeneracy point 
(6 = it/2) is always rounded out by higher order terms, 
and the anisotropy energy is an even analytic function 
of cos(#). The degree of rounding is a complex many- 
body issue, which in our DFT calculations is affected 
by the smearing of occupation numbers at the HOMO. 
To go beyond our calculations, one would have to em- 
ploy a true many-body framework, where the absence of 
the otherwise necessary smearing parameter would imply 
even larger anisotropics. 

From fig. [5] it is clear that the total Chern number of 
4 originates in equal parts from spin and orbital angular 
momentum. Without SOI there is no angular momen- 
tum contribution and in this case the total Chern num- 
ber is 2 coming from the spins in the (i-shells. With SOI 
the quasi-degeneracy is lifted along the easy (symmetry) 
axis and the total energy is lowered via the mixing of d xy 
and d x 2_ y 2 into axial eigenstates with mi — ±2 maxi- 
mal splitting and a continuous electron path in a plane 
perpendicular to the easy axis magnetization direction. 
At the equator, the 5* levels approach each other and 
the curvature becomes quasi-singular. Here the d xy and 
d x 2_ y 2 characters separate and we find that the HOMO 
in the a;-direction (<j) = 0) is of pure d x 2_ y 2 character (see 
fig. 0(c) and (e)). 

The effect of the equatorial curvature ridge in our vari- 
able transformation of the effective big-spin Hamiltonian, 
is to alter the local metric (|12[) and repel grid-points from 
the equator in order to increase the subtended curvature 
areas in this region. This leads to a transformed mag- 
netic anisotropy landscape that has a much wider bar- 
rier (see fig- El) • Such an enormous increase in effective 
barrier width should definitely lead to experimentally ob- 
servable effects in terms of spin-lifetimes and tunneling 
probability. 

The Chern number is 4, which tells us that we may use 
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FIG. 9: (Color online) Eigenlevel structure around the 
HOMO and atomic site projected orbital character, for the 
easy (x for Ni2 and Pd2, z for Pd^) and the hard (z for N12 and 
Pd2, x for Pd^) direction, as well as over the unit sphere for 
the HOMO. The orbitals in parenthesis are the same as the 
ones displayed, only with an equatorial ^-dependence that is 
in relative anti-phase. 

spherical harmonics of up to order L = 8 in our quanti- 
zation process. It is interesting to note that we actually 
need all of these to accurately describe the transformed 
magnetic anisotropy landscape. Omitting the smallest 
8:th order coefficient leads to the formation of a non- 
physical shallow local minimum along the equator. 

Paralleling the close relationship between Co and Rh, 
Pd and Ni have similar valence electron configurations 
- the 4d elements behave like their 3d counterparts. In 
the absence of SOI, both C02 and Rli2 have a molecular 
groundstate A(S — 2) s and a Chern number J = 4 with 
SOI. Pd2 and Nig both have a Chern number of 1 (with 
and without SOI) with a S(5 = l) s ground state. Their 
anisotropy landscapes are quasi-easy planes perpendicu- 
lar to the dimer line, although the anisotropy energies 
are small - balancing rather delicately between easy axis 
and easy plane. The curvature is essentially constant 
and equal to the Chern number, wherefore no topolog- 
ical effect enters via the variable transform (see fig. [7|). 
The combined s-shell in the dimer, demotes two atomic 
s-electrons to the d-shell, leaving only two d-holes in the 
dimer. The net d-spin contribution yields the Chern 
number 1, with no orbital contribution. Pd2 has a much 



larger dimer separation than M2, which results in a dif- 
ferent ordering of the levels. Fig. [5] (d) reveals that the 
7r*-level is lower in energy than the 4s(cr g ) in Pd2, result- 
ing in a 7T* HOMO, whereas Ni 2 retains the 5* (fig.[S](a)). 
In both cases, the doubly degenerate HOMOs are dou- 
bly occupied. The antipodal shifts now cancel out just 
like for most of the sub-levels. This highlights the sit- 
uation for which high anisotropy and strong topological 
effects are present: when there is a doubly degenerate and 
singly occupied HOMO in the absence of SOI. The quasi- 
degenerate occupied levels in Pd2 and N12 both have a 
highly singular level curvature, but these also cancel out, 
resulting in a trivial total curvature with very small vari- 
ations around the Chern number. It is clear that the 
case of a singly occupied, doubly degenerate HOMO is 
the most natural candidate for a system in which there 
is a powerful topological effect. However, Pd^ points to 
another possibility. In reducing the dimer distance - the 
s-level drops in energy, until it actually intersects the 
HOMO-LUMO gap. In doing so, the strong topological 
effect (see fig. |SJ) is triggered due to the mixing with the 
7r* doublet (see fig. 0(g)). The non-trivial curvature en- 
ters the transform and causes the effective barrier for Pd^ 
to widen and the effective Hamiltonian acquires a fourth 
order contribution (see table UV[). Note the inverse re- 
lation between the s-character of the HOMO (fig. O (h)) 
and the curvature (fig-EJ), i.e. the curvature barrier grows 
as the s-character diminishes and the d a;z ((ij /2 )-character 
increases. 

The situation observed in Pd^, where an accidental 
degeneracy triggers a topological effect, by intermixing 
with the HOMO-LUMO doublet, can also occur in re- 
laxed systems. For example, in C^r^, there is the ex- 
act same mechanism as in Pd^ present at the HOMO, 
where an s-level causes an accidental degeneracy with a 
7r HOMO, for a wide range of dimer distances. 



F. C03 and Co| 

Just like the symmetrized C045, the equilateral trimer 
C03 has a highly degenerate level structure and are there- 
fore unstable against Jahn- Teller distortions (see fig-d]). 
Fig. [TU] (a) shows the levels in close vicinity of the HOMO 
in the hard (x) and easy (z) direction for C03 and the cor- 
responding directions for the Jahn- Teller-distorted C03, 
where the hard and easy direction is interchanged relative 
the equilateral structure. 

For the C03 the orbital mixing changes very little be- 
tween the two directions and there is a large HOMO- 
LUMO gap of 0.03-0.04 eV. By contrast, the Co^ has a 
quasi-degeneracy in the hard direction of approximately 
4 meV that is lifted to around 35 meV in the easy di- 
rection. The quasi-degeneracy that is lifted by the Jahn- 
Teller distortion, induces large variations from the Chern 
number in the curvature (see fig. [5] and [7]). In fig. [3] we 
see that the level Chern number for the HOMO is -3/2, 
meaning that the HOMO level curvature that goes quasi- 
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FIG. 10: (Color online) Voxel projections of the C03" HOMO 
and LUMO charge densities for different directions of the 
magnetization, (a) shows the energies of the HOMO and 
the LUMO for the isosceles and the equilateral trimer in the 
hard and easy direction, (b) shows the HOMO charge density 
for the equilateral trimer in a top view using maximal inten- 
sity projection. Here, a small amount of d-charge forms axial 
states and the HOMO-LUMO splitting is maximal, lowering 
the total energy, (c) and (e) show the HOMO charge density 
when the magnetization is in the hard plane along the equa- 
tor. Here, the d-character separates, occupying the same sites 
as the p-orbitals. (d) and (f) show the corresponding LUMO 
charge densities which occupy increasingly disjoint regions of 
the cluster as the gap closes in. 
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FIG. 11: (Color online) The level structure around the HOMO 
for C05 and Co3Fe2 and voxel projections of the Co3Fe2 
HOMO charge density for different magnetization directions, 
(a) shows that the substitution of the axial Co atoms in C05 
for Fe closes the HOMO-LUMO gap from approximately 0.4 
to 0.002-0.01 eV, causing an extreme quasi-degeneracy and 
impact on the Berry curvature, (b) shows a maximal inten- 
sity projection of the HOMO charge density when the mag- 
netization direction is parallel to the bi-pyramidal easy-axis 
(top view), (c) and (d) show the HOMO densities with the 
magnetization in hard plane. The variations in d-character 
are responsible for the equatorial oscillations in the curva- 
ture. In (e) and (f) the LUMO charge densities when the 
magnetization is in the z- and x-direction are shown. 



G. Co 5 and Co3Fe 2 



singular has negative curvature. Although this negative 
curvature is shifted by the accumulated Chern number 
coming from the occupied sub-levels (level 26 has an ac- 
cumulated Chern number of 3 in fig. EJ), negative curva- 
ture residues remain, implying that we can only approx- 
imately describe the topological impact by omitting the 
negative dips in the transform. 

Fig. [TU] shows voxel projections of the HOMO and 
LUMO charge densities for different directions of the 
magnetization. When the magnetization is in the easy 
(z) direction, perpendicular to the trimer plane, each 
atom forms a different combination of on-site d- and p- 
character, symmetrically around the z-axis (see fig. [10] 
(b)). Here, a small mixture of approximately 10% d xy 
and d x 2_ y 2 in the whole system is responsible for the 
lowering of the HOMO energy. On the equator (see 
fig. [10]) (b)-(f)), these separate, different atoms are oc- 
cupied by the HOMO charge and the HOMO-LUMO be- 
come quasi-degenerate. The splitting increases slightly 
along the equator as the magnetization changes to <fi — ^ 
(see fig. [10] (e)), enough to produce the observed oscilla- 
tions in the curvature. In the hard direction, the LUMO 
density occupies increasingly disjoint regions of the clus- 
ter, as the HOMO-LUMO gap closes in. The occupied 
configuration where the atoms are unequally occupied 
is higher in energy than the symmetric configuration in 
along the z-axis, resulting in a bistable configuration with 
a total magnetic anisotropy energy of 3 meV. 



Diagrams of the level structure around the HOMO for 
Co3Fe2 and C05 are shown in fig. [TTJ(a). When the axial 
Co atoms in C05 are substituted for Fe, the character of 
the HOMO changes from p- to d-like. In addition, the 
HOMO-LUMO gap decreases from approximately 0.4 eV 
in C05 to 2-10 meV in Co3Fe2. This dense level structure 
in Co3Fe2 causes a level-crossing when SOI is turned on, 
changing the Chern number from 13/2 to 11/2. In C05 
however, there is no level-crossing, and the Chern number 
does not change with the inclusion of SOI. 

Fig.[5]shows that the substitutional Fe atoms compress 
the cluster along the bi-pyramidal axis and the center 
triangle increases in size. This points to the presence of 
an attractive force between the Fe atoms that is coun- 
teracted by the Co triangle, causing a more symmetric 
system with an accidental degeneracy at the HOMO. The 
anisotropy/atom (see table ITT|) increases by a factor five 
with the substitution - largely due to the uncompensated 
HOMO-downshift. Pederson et. Al^£ report a match- 
ing anisotropy/atom of 0.1 meV and moment of 13 [Ib 
for C05 with an increase in anisotropy/atom with a fac- 
tor of approximately 5 in Co3Fe2- The moment of 13 
Hb matches our non spin-orbit Chern number, but this 
changes as SOI is included (see table [TJ. Since the large 
fluctuations in the curvature at the equator (see fig. [5]) 
cancel out to a large extent in the 0-variable transform 
(fTTj) . and since there is no ^-dependence in the anisotropy 
landscape, the prevailing topological effect can be cap- 
tured in the transform by cutting the negative dips (see 
fig.©. 

Fig. [TJJ (b)-(d) show maximal intensity projections of 
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FIG. 12: (Color online) Energy levels and charge densities for 
the HOMO and the LUMO in Mn 2 N 3 . (a) shows the HOMO 
and LUMO energies in the hard (z) and easy (x) direction. In 

(a) , (b) and (c), the total charge density, the magnetization 
density in the FM variation and in the AFM variation are 
shown, respectively, (e) and (f) show the charge density of 
the HOMO and the LUMO in the ^-direction and (g) and (h) 
in the a;-direction. 

the HOMO charge density along the easy (z) axis and 
along two symmetry directions along the hard plane (the 
Co3-plane). These images reveal that there is not much 
HOMO charge on the Fe atoms and that the dynamics 
of the HOMO is governed by the Co3-triangle. In Fig. [TT1 

(b) we see that the Co do not really form complete axial 
states, rather the high anisotropy comes from the fact 
that the states corresponding to the hard plane (Fig. [11] 

(c) and (d)) are relatively higher in energy. The variation 
of the d-orbitals on the different Co when the magnetiza- 
tion changes in the hard plane is responsible for the large 
equatorial oscillations. In fig. QT] (e) and (f) the LUMO 
charge densities are shown for a magnetization in the z 
and cc-direction respectively. The z LUMO charge den- 
sity has maximum on the Fe atoms, and is much higher 
in energy than the symmetric Co triangle in the HOMO. 



H. MrizNy 

Finally, we turn our attention to the hybrid Mn^Ny 
clusters, partly motivated by the work of Hirjibehedin 
et. Al3i. In order to properly simulate their system, it 
is necessary to include some layers of Cu, which, although 
beyond the scope of this article, certainly warrants fur- 
ther investigation. As the interesting topological effects 
occur only for the AFM configurations, we shall focus 
primarily on those. 

We start with the Mn 2 N 3 -clusters. In Mn 2 N^ the 
atoms have been locked in a configuration that corre- 
sponds to their positions as they are embedded in the 



covalent network of the surface. This cluster does not 
differ much from the Mn 2 N3, where the Mn have been 
frozen at a distance of 3.6 A and the N where allowed to 
relax starting from off Mn-axis positions. The angle of 
the endpoint N is some what steeper and the distance to 
the Mn longer because of the presence of the surface Cu 
atoms. 

Fig. [H(a) show how the HOMO-LUMO levels vary be- 
tween the easy (x) and the hard (z) direction, for M112N3. 
In Mn 2 N3 the gap varies between 7 and 0.5 meV between 
the hard and easy direction, where as the Mn 2 Ng has the 
same qualitative shifts but with a smaller variation from 
3.5 to 1 meV. In Fig. [H (b) the total charge density is 
shown, where the interspacing N atoms hold much less 
valence charge, (c) and (d) show the magnetization den- 
sity in the Mn-axis direction, for the FM and the AFM 
variation respectively. The fact hat the N atoms magne- 
tize oppositely the Mn atoms in the FM configuration, 
is an indication that the N stabilize the AFM configu- 
ration. The FM and the AFM configurations are rather 
close in energy (see table [H]) although the AFM is lower 
in energy by the order of ten meV. As noted in by Rao 
et Al4&, the presence of N in small Mn clusters, dramati- 
cally increases the binding energy and generally enhances 
the Mn magnetic moments, highlighting the role of the 
N as mediators of the AFM coupling. We find a Chern 
number of 1/2 in the AFM and 5/2 in the FM variation, 
which indicates that the N contribute to the effective 
cluster spin dimension. This configuration could change 
by the inclusion of surrounding surface Cu atoms. In the 
case of the dimer, a Heisenberg Hamiltonian fitted to ex- 
perimental data does not lead to distinguishable results 
between different effective atomic Mn spins. 

Fig. [H ( e ) and (g) show the HOMO and (f) and (h) 
the LUMO charge density, in the hard (z) and easy (x) 
direction respectively for the Mn 2 N3 (Mn 2 Ng follows a 
similar pattern, although less pronounced). The p-like 
charge on the N connects with the d-like Mn charge in one 
end of the cluster. The charge configurations in fig. [H 
(e) and (f) results in a higher energy for the LUMO 
(level Chern number —1/2) and a lower energy for the 
HOMO (level Chern number 1/2). As the magnetiza- 
tion is changed to the x-direction, the quasi-degenerate 
HOMO and LUMO occupy completely disjoint regions 
of the cluster. The energy difference between these sym- 
metric mirror configurations is very small, resulting in 
severe quasi-degeneracies. These occur perpendicularly 
to the cluster plane and appear in the form of Gaussian 
shaped quasi-singularities in the curvature (see fig. [|5]). 
Their effect in the transform is to compress the low bar- 
riers to a thin line. After quantization this yields an ef- 
fective quasi-plane Hamiltonian with reduced anisotropy 
(see tab. [TV]). The HOMO shift is not large enough to 
turn the z-direction into the easy direction, it does how- 
ever reduce the anisotropy somewhat. In this case the 
MAE originates in the lower levels, whereas the its non- 
trivial topology is coming from the HOMO-LUMO level. 
The quasi-degeneracy that appears is the most severe one 
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FIG. 13: (Color online) Densities for MnaN^. (a) shows the 
total charge density and (b) the magnetization density in the 
easy (z) direction, (c) and (d) show the HOMO charge den- 
sities in the easy and hard direction, with perspective views 
in (e) and (g). (f) and (h) show the LUMO for the relaxed 
cluster which is 0.3 eV above the HOMO. 

that we have encountered (see fig. [5]). In view of this, 
here we expect that our single spin procedure to break 
down. Nevertheless, the transform gives a good indica- 
tion of what is going on - the anisotropy is basically be- 
ing erased by the topology of the system. The absence of 
spin-flip excitations in the conductance spectrum^ could 
very well be a manifestation of this powerful topological 
effect, even when the Chern number is different from zero. 

We now address the Mn 3 N2. The AFM variation with 
Chern number 5/2 has a much lower energy than the 
FM variation with Chern number 15/2. This is consis- 
tent with representing the Mn with atomic spins of 5/2 
and the anisotropy is of the order of what one would ex- 
pect from the conduction spectrum. M^N^ sticks out 
as there is no obvious correlation between the inverse 
squared of the HOMO-LUMO gap and the Berry curva- 
ture (see fig. [5]) . The gap size does not vary that much 
between hard and easy directions and is quite large, be- 
tween 20-40 meV. In addition, only the HOMO moves 
as the magnetization direction is changed. The HOMO 
up-shift is in fact compensated by the sub-HOMO level, 
whereas the LUMO does not move at all. This is the 
reason for the small anisotropy the odd appearance of 
the curvature. Fig. [13] (a) shows the charge density and 
(b) the magnetization density in the easy-direction (along 
the cluster axis), revealing no magnetic charge on the N 
atoms, (c) and (d) show the HOMO charge density in 
the easy- and hard-direction respectively. Varying 6 from 
to 7T, we see that the Mn p- and d-character traces out 
a contour like the one observed in the curvature. The 
curvature causes a slight widening of the blocking bar- 
rier (see fig. [7]) with a small fourth order contribution 
in the effective Hamiltonian (see tab. IIV[) When M^N^ 
is allowed to relax into Mn 3 N2 (see fig. H]), the N po- 
sition themselves relatively closer to the endpoint Mn, 
and the HOMO now shows no charge on the center Mn. 
The endpoint Mn and the N now form orbital configura- 
tions where the LUMOs are approximately 0.3 eV above 
the HOMO, resulting in a trivial and constant curvature. 



This highlights the crucial role of the epitaxial registra- 
tion and the symmetry enforced by the surface. Adding 
two N atoms on the endpoint to create M^N^ we find 
that the Chern numbers change to 9/2 for the FM and 
3/2 for the AFM variation, consistent with instead repre- 
senting the Mn with spin 3/2. Clearly the extra N alter 
the Mn electronic configuration, pointing to a sensitive 
dependence on the endpoint N positions. The FM vari- 
ation yields a quasi-easy plane with a lower global min- 
imum, and only in the AFM variation of Mn 3 N;![ do we 
find an easy axis parallel to the cluster axis. Relaxing the 
atoms on a line again returns the AFM to the quasi-easy 
plane configuration. The Mn 3 N| easy HOMO charge ex- 
hibits ^-character on the endpoint N, but the hard di- 
rection HOMO charge has a dominant d^y-character on 
the center Mn that is 17 meV higher in energy. This level 
is very close in energy to the hard direction LUMO, with 
a d 2 .2_ 2/ 2-character - just a rotated version of the hard 
HOMO. The non-trivial Berry curvature manifests itself 
precisely because of this rotational symmetry. The large 
shift in the HOMO of 17 meV is enough to tip the scales 
and change the MAE from quasi-easy plane in Mn 3 N4 
to an easy axis in Mn 3 N^. The quasi-degeneracy in the 
hard plane causes a widening of the blocking barrier (see 
fig. [7]). Quantizing with J = 3/2 entails an averaging into 
the allowed spin-space, such that the widened barrier is 
represented by a greater effective barrier height. 



IV. ALTERNATIVE APPROACH AT A 
DEGENERACY POINT 

The extensive numerical studies presented above show 
that close to a HOMO-LUMO degeneracy we face the 
most interesting and difficult situation, where our giant- 
spin formalism is pushed to its limit: on one hand the 
(partly avoided) level-crossing signals that Berry cur- 
vature effects are important and large Berry curvature 
fluctuations modify the effective energy anisotropy lans- 
dcape. On the other hand, when the HOMO-LUMO gap 
becomes too small, the ensuing singularities in the Berry 
curvature render the implementation of our formalism 
problematic, particularly when the total curvature be- 
comes negative and the variable transformations given 
in Eqs. [Til EH are not valid. In this case the descrip- 
tion of the low-energy dynamics in terms of one single 
"giant-spin" breaks down, and it is necessary to reintro- 
duce explicitly the relevant electronic degrees of freedom 
responsible for the non-adiabaticity. We sketch here an 
alternative treatment for the spin dynamics that should 
remain valid when a true electronic degeneracy is present. 
Our idea is based on an analogy of our problem with the 
Jahn- Teller effect (JT) in the vibronic studies of poly- 
atomic molecules. (For a review of the JT effect see 
Ref. H^.) In the static JT effect, an electronic level de- 
generacy, usually present for a very symmetric molecular 
configuration, is lifted by a distortion of the molecule 
to a lower-symmetry configuration, accomplished via a 
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coupling of the electronic variables with the nuclear de- 
grees of freedom. The molecular symmetry, reduced in 
the static JT effect, is restored by coherent tunneling 
between equivalent distortions, a phenomenon known as 
dynamical JT (DJT) effect. 

The coupled dynamics of the (fast) electronic degrees 
of freedom and the (slow) nuclear motion is usually 
investigated in two steps within the adiabatic Born- 
Oppenhcimer (BO) approximation. Here the electronic 
problem is first solved for a fixed arbitrary configuration 
of the nuclear variables, assumed in the first instance to 
be classical. The resulting electronic energies act as po- 
tential energy functions for the nuclear dynamics prob- 
lem, which is then treated quantum mechanical. It turns 
out that near an electronic degeneracy giving rise to the 
JT phenomena, non-adiabatic corrections to this proce- 
dure in the form of Berry phases must be taken into ac- 
count. The important role of such non-adiabatic correc- 
tions was already recognized in early studies of electron- 
lattice interactions in molecules (see Ref. for a re- 
view of this work) . However, it was after Berry's seminal 
work22 that these corrections were recognized to be ex- 
amples of the ubiquitous quantum geometric phases of- 
ten arising in scmiclassical treatments of complex quan- 
tum systems containing slow and fast degrees of freedom. 
The study of DJT effects in molecular systems remains a 
topic of great interest, which presently includes investi- 
gations of novel nanostructures carried out with the help 
of powerful electronic structure calculations. In this con- 
text it has recently been pointed oui£i that electronic 
structure calculations generating global minimum con- 
figurations and associated vibrational frequencies might 
be faulty if Berry phase effects, when present, are ne- 
glected. In general Berry phase terms affect profoundly 
the quantum spectrum of the vibronic dynamics. 

In order to construct our analogy with the DJT ef- 
fect, we first review the classical case of the JT E <g> e 
problem (a doubly-degenerate electronic E term interact- 
ing with doubly-degenerate e vibrational mode) 25 i 52 ' 53 . 
The Hamiltonian for the two electronic states {|a),|6)} 
and the doubly degenerate normal vibrational coordi- 
nates Qi = pcosd, Q2 — psinO is H = Hq + Hi, 
where H = ^2 i=ajb fedho with h 2 &h {p,8) being 

the Hamiltonian of a two-dimensional isotropic harmonic 
oscillator describing the free normal modes. Hi — 

(k pe~ t6 + \g p 2 'e l26 ^j |i)(i|+H.c. is the vibronic interaction 
Hamiltonian. The coordinates Qi=Q2 = 0^p = 
represent the symmetric (undistorted) configuration of 
the molecule. The adiabatic BO potential surfaces E± 
are obtained by diagonalizing the electronic part of the 
H; that is, the nuclear kinetic energy is set to zero and the 
vibrational coordinates are treated as classical variables. 
For the linear E <£> e JT problem, for which g — 0, E± — 
^p 2 ikp, and the corresponding adiabatic electronic wave 
functions are |±(0)) = (e- %9 / 2 \a) ±e l8 / 2 \b))V2. The low- 
est potential energy surface has a continuous minimum at 
po = k, which gives rise to the JT effect. The electronic 



wave functions | ± (£?)) are double- valued when trans- 
ported adiabatically around the degeneracy in the elec- 
tronic spectrum {E+ = EJ) occurring at p = 0. The sign 
change in |±(0)) upon encircling the origin is a manifesta- 
tion of the topological Berry phase arising from the con- 
ical intersection of the adiabatic potential surface. The 
sign change has dramatic effect on the quantum vibration 
spectrum. The simplest way to see this is to consider the 
large k regime, where the two potential surfaces are well 
separated. We can therefore quantize the nuclear motion 
by restricting to the lowest energy surface E-. The rel- 
evant potential energy to be added to /i2dho is therefore 
V = E_ + H BH , where H BH = \ £4=1,2 9{ ^ )l %^ is 
the Born-Huang term 2 ^^^^. 

Since the total vibronic wave function ty(p,6)- = 
\p,6)~ <8 I — (#)) must be single- valued, the nuclear part 
\p,0)- must change sign when going around the origin 
to compensate the double- valuedness of | — (8)). This 
boundary condition causes the spectrum of /i2dho to have 
half-odd integral vibronic angular momentum quantum 
numbers. Alternatively, the sign change of \—){8) can 
be absorbed by re-defining its phase. This phase change 
introduces a fictitious magnetic vector potential in the ki- 
netic energy of the nuclear motion, known as Mead-Berry 
vector potential or Mead-Berry connection, whose curl is 
the Berry curvature. The presence of this fictitious mag- 
netic field causes the same half-integer quantum numbers 
for the vibronic spectrum. 

Quadratic and higher coupling terms in Hi are in gen- 
eral important and give rise to more complex adiabatic 
potential energy surfaces. In particular, other conical in- 
tersections are possible, apart from the one at p = 0. 
Typically in the E ® e problem three other degeneracies 
are possible, each of them generating a Berry-Mead vec- 
tor potential. In this case the overall behavior of the elec- 
tronic wave functions and the resulting vibronic spectrum 
depend crucially on the dominant tunneling paths con- 
necting potential energy minima, around these degenera- 
cies: if paths encircling all four degeneracies are possible, 
the overall Berry phase is 47r, that is, the wave functions 
are single- valued and the vibronic quantum numbers are 
integers. 

In our problem with magnetic clusters, the slow nu- 
clear motion is replaced by the collective spin degree of 
freedom of the cluster S — Sh. This is couple to two or 
more electronic levels, including the HL levels, which are 
degenerate for given spin orientations. It is important 
to clarify that in this case we are not interested in the 
nuclear dynamics: the molecule configuration is assumed 
to be fixed. The Hamiltonian for the collective spin and 
the HL electronic levels {\H), \L)} is 

H = H + Hi , (22) 

where 

H *= E (^ + U 0))\i)(i\, (23) 

i=H,L 
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with V(S) being the anisotropy energy functional with- 
out the contribution of the HL levels. Here we take 5 
to be the total spin moment of the cluster when SO cou- 
pling is switched off. The term Hi describes the coupling 
between the collective spin and the HL electronic levels, 
and is written as 

H 1 = £ V t (S)\ t )(i\ + (v HL (S)\H)(L\+H.c) . (24) 

i=H,L 

The functions Vl, Vh, and Vrl = V£ H have to be cho- 
sen so that the eigenvalues E± (5) of Hi reproduce the 
HL energy landscapes, with the associated intersections 
as a function of the orientation of 5. The eigenvalue and 
eigenvectors of Eq. [22] can be obtained numerically by 
expanding in the basis set | ± 5) <8> \m,S, where |m, S 
are eigenvectors of 5 2 and S z . However, since we want 
to single out the effect of the Berry phase on the slow 
variable dynamics, it is useful to think in terms of the 
adiabatic semiclassical approach and use an approximate 
spin Hamiltonian that includes the lowest energy surface 
surface £L (5) only. Like for the D JT problem, a conical 
intersection will generate a Berry phase in the electronic 
wavefunction, which in turn will affect the quantum spec- 
trum of the collective spin variable 5. In the presence 
of several conical intersections, the overall effect of the 
Berry phase on the electronic wave function (and there- 
fore on the dynamics of the collective spin) will depend 
on the character of the dominant tunneling paths con- 
necting nearby energy minima through different saddle 
points around the degeneracies. As an example we can 
consider the case of the Co trimer. When the cluster 
is fixed in the static JT distorted configuration, there is 
no intersection between the HL energy function. There- 
fore the Berry phase is absent and the total spin 5 is 
the same as for the case in which SO is absent. However 
when the trimer is fixed in a symmetric configuration, the 
HL gap vanishes for a few spin orientations in the plane 
of the clusters. Note that, as in the ordinary static JT 
effect, the system lifts the degeneracy by choosing a "dis- 
torted" configuration of the slow variables (away from the 
"symmetric" one) that lowers the total energy. In C03 
this configuration is the spin orientation orthogonal to 
the trimer plane. An analysis of the HL energy surfaces 
in the trimer plane reveals that there are 6 degeneracy 
points separated by very shallow energy barriers. The 
fact that the Chern number of the cluster changes from 
5 = 5/2 to 5 = 3/2 when SO is included, seems to indi- 
cate that paths that encircle an even number of conical 
intersections have dominant tunneling rates. 

V. CONCLUSIONS AND OUTLOOK 

The crossover between macroscopic magnetic particles 
and atomic spins passes through the rich and complex 
world of magnetic nano-clusters. Interest in exploring 
this world is growing because of the push toward denser 



magnetic information storage media and an ensuing in- 
terest in achieving a deeper understanding of fundamen- 
tal limits, and in identifying potentially attractive sys- 
tems. From the classical point of view the most im- 
portant property of a magnetic particle is the depen- 
dence of its energy on moment orientation, i.e. its mag- 
netic anisotropy. Magnetic anisotropy in magnetic nano- 
clusters is due almost entirely to spin-orbit interactions 
which induce an orbital contribution to the moment and 
therefore a dependence of energy on the orientation of 
the moment relative to the spatial arrangement of atoms 
in the cluster. 

In this paper our interest has been focused on the 
quantum effects which become important in small mag- 
netic clusters. In particular, we have explored a SDFT- 
based approach to derive a giant spin low-energy effective 
Hamiltonian from microscopic physics. This Hamiltonian 
is a quantum generalization of the magnetic anisotropy 
energy, and describes a well-defined group of low-lying 
energy levels which can be identified with the magnetic 
moment orientation degree-of-freedom. 

Our approach is orbital-based, rather than local- 
moment based, and is intended to enable a practical de- 
scription of systems in which the spins are carried by 
orbitals which are widely spread over many of the clus- 
ter atoms. Loosely speaking, our approach is intended to 
enable a quantum description of metallic magnetic nano- 
clusters in which the spin-moment per atom in the ab- 
sence of spin-orbit interactions is not an integer. The 
quantization rules of the classical anisotropy energy are 
calculated by adding the adiabatic Berry phase contribu- 
tions from all itinerant orbitals which behave collectively 
in the giant spin approximation which is appropriate for 
slow (low-excitation energy) dynamics. We believe that 
the approach we have taken will be accurate whenever 
the fundamental assumptions which underly the giant 
spin approximation are valid. 

In our approach both the anisotropy energy functional 
and the Berry curvature functional are calculated by solv- 
ing the DFT Kohn-Sham equations using standard spin- 
polarized density functional methods that include the 
spin-orbit interaction. The average of the Berry curva- 
ture over magnetization directions, a topological invari- 
ant known as Berry-phase Chern number, determines the 
total angular momentum J and reduces to the spin an- 
gular momentum 5 in the absence of spin-orbit interac- 
tions. In our theory the Chern number determines the 
dimension of the quantum giant-spin Hamiltonian. We 
find that J can deviate considerably from the total spin in 
the absence of spin-orbit interactions 5 when the HOMO- 
LUMO gap is small. This is often the case in transition 
metal magnetic nano-clusters, which are often charac- 
terized by a dense level structure near the Fermi level 
even when the clusters are quite small. When the HL 
gap is small, the magnetization-dependent Berry curva- 
ture has large fluctuations and its effect on the effective 
spin Hamiltonian, which is usually neglected in treat- 
ments of molecular magnets, is very significant. We call 
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this effect topological since it is related to the topolog- 
ical Berry-phase Chern number, and can be intuitively 
described as a stretching or a compression of the clas- 
sical magnetic anisotropy landscape. Although we have 
not performed calculations for these systems, we do not 
expect the topological effect to be important in molecu- 
lar magnets with well defined local moments. For metals 
however, the Berry phase effect modifies qualitatively the 
collective dynamics of the magnetization degree of free- 
dom. It is clear that extremely large topological effects 
often signal a breakdown of the giant-spin approxima- 
tion. In that case we view our approach as an attempt 
to construct the best possible giant-spin model. 

We have carried out our DFT-based quantum giant- 
spin model construction procedure for dimers and for 
other clusters containing up to five identical atoms. We 
find that dimers can display a huge magnetic anisotropy 
because of their strong axial symmetry. For example we 
find anisotropies as large as « 500 K in Rh 2 , which are 
accompanied by important topological effects that widen 
the magnetic anisotropy barrier between the two equiva- 
lent uniaxial minima. Both properties are closely related 
to the HOMO-LUMO degeneracy mentioned above, and 
suggest that magnetic dimer systems could be of poten- 
tial technological importance for magnetic storage ap- 
plications if they could be embedded in an appropriate 
lattice. Such large magnetic anisotropies for such small 
systems could represent the ultimate limit of magnetic 
storage capability. Furthermore thicker energy barriers 
will tend to increase the spin-relaxation time and sup- 
press the microscopic quantum mechanical tunneling of 
the magnetization. Larger clusters such as C03, C04, and 
C05 , typically manage to avoid quasi-degeneracy between 
the HOMO-LUMO levels by means of Jahn-Teller distor- 
tions. In this case the impact of the topological Berry 
curvature of the spin Hamiltonian is minor. If however a 
degeneracy if forced on the system by imposing a particu- 
lar symmetric configuration, the singular Berry curvature 
once again plays a dominant role. Here we have shown 



that the Berry curvature can even become negative and 
our procedure, strictly speaking, breaks down, signaling 
that a description of the low-energy spectrum in terms 
of an individual collective spin is not possible. For these 
cases one is compelled to reintroduce the relevant elec- 
tronic degrees of freedom involved in the degeneracy at 
the Fermi level. We have sketched how this procedure 
should be carried out in Sec. |IV] 

We have also looked at several hybrid clusters, in par- 
ticular at linear Mn^N^-clusters in which the interatomic 
distances have been locked to mimic the effect of the epi- 
taxial registration on a CuN surface. The study of this 
system is partly motivated by recent STM experiments 20 
that are able for the first time to probe the collective spin 
dynamics of quantum engineered magnetic clusters via 
inelastic electron tunneling. Experimentally these clus- 
ters are found to order antiferromagnetically. We have 
investigated this possibility and carried out our proce- 
dure for this situation. We find that these systems show 
different topological effects depending on the configura- 
tion studied, and exhibit high sensitivity to the position 
of the N atoms coming from the CuN substrate. The re- 
sults offer an alternative interpretation for the Mn dimer 
differential conduction spectrum and in the case of the 
Mn trimer our findings are compatible with the exper- 
imental results. We should however point out that our 
calculations are for free clusters only, namely we did not 
include any complicating effects originating from the sup- 
porting CuN surface that are likely to be important. 

Our findings indicate that powerful effects of non- 
trivial spin-space topology are present in many small TM 
clusters, and cannot be ignored when determining the 
collective spin dynamics. Furthermore, our approach of- 
fers a way to extract the total spin of the system, without 
assigning locally defined atomic spins. Only more de- 
tailed calculations that include supporting surfaces and 
comparison with more accurate experiments will tell if 
the effects predicted by our theory can be observed in 
these systems. 
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